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Abstract 

A recently introduced numerical approach to quantum systems is analyzed. The basis of a Fock 
space is restricted and represented in an algebraic program. Convergence with increasing size of basis 
is proved and the difference between discrete and continuous spectrum is stressed. In particular a 
new scaling low for nonlocalized states is obtained. Exact solutions for several cases as well as general 
properties of the method are given. 

1 Introduction 

Recently an attractive possibility of modelling M-theory through relatively simple quantum me- 
chanical systems [1] has occurred. They emerge from the dimensional reduction of supersymmetric 
gauge theories and provide a simple laboratory to study many properties of supersymmetry [2,3] . 
It follows from [1] that there is an strong connection between M-theory and SU{Nc — > oo) super- 
symetric Yang-Mills quantum mechanics (SYMQM). However supersymmetric quantum mechanics 
have much longer history. Various schemes have been analyzed to try to solve the hierarchy problem 
including the idea of breaking SUSY. This was the reason why SUSY was first studied in the simplest 
case of quantum mechanics (SUSYQM)[2]. Apart from its physical meaning SUSYQM gave also a 
deeper understanding of why certain potentials are analytically solvable and other are not [4]. The 
SYMQM gauge systems were studied for the first time in [3] where the exact spectrum including the 
ground state of SYMQM D=2 was given. Later on the extension for arbitrary SU(N) gauge group 
was also obtained [5]. SUSYQM is known to have continuous spectrum due to the fermion-boson 
cancellation [6]. According to BFSS hypothesis there should be a bound state at the threshold of 
the spectrum. However, since there are no exact solutions one is forced to use numerical methods. 



In this paper we discuss in details a numerical approach of solving quantum mechanical systems 
proposed in [7,8] and already investigated in [9-13[. Next section contains formulation of the method 
as well as its general properties. We introduce a cutoff, N, and by means of an algebraic program 
analyze a complete dependence of the spectrum on the cutoff. We prove that the eigenvalues converge 
towards exact (i.e in the infinite Hilbert space) spectrum. In section 3 we give the exact spectrum 
of the momentum and coordinate operators at arbitrary finite A'^. The asymptotic behavior with 
A'^ — > CO is derived in section 4 where a new scaling law, required to recover the infinite Hilbert space 
limit, is formulated. The scaling and its universality is discussed in section 5 by giving the exact 
spectrum of a free particle in quantum mechanics. Interestingly, this solution differs only a little in 
comparison with the eigenvalues of the hamiltonian for D — 2 supersymmetric Yang-Mills quantum 
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mechanics at finite cutoff [14]. We prove that the continuum spectrum in quantum mechanics gives 
rise to the power-like dependance on the cutoff. This result is important in studying supersymmetric 
systems where the distinction between continuum and discrete spectra is an important issue. In 
section 6 we use numerical data in order to verify the theoretical results. The implementation of the 
approach in Mathematica code will be discussed there in details. 



2 A cut Fock space 

Every quantum hamiltonian can be represented in the eigenbasis of a harmonic oscillator 

{\n) = ^\0),nGN}, (1) 

where a, are the normalized annihilation and creation operators respectively. The correspondence 
between a, and Q, P (coordinate and momentum operators respectively) reads 

Q=^^ia + a^P=^ia-a^ (2) 

Since this basis is countable it is very convenient to use it in numerical applications. One can limit 
(1), e.g. n ^ N, then calculate the finite matrix representation of any hamiltonian and numerically 
diagonalize above finite matrix to obtain complete spectrum and the eigenstates of the system ^. 
The procedure is simple and essentially numerical, however a number of theoretical questions arises 
while analyzing it. They will be discussed in this paper. 



We denote 

H^'^^ = [H]ij = hij i,j = l,...,N + las operator H in a. cut Fock space (cutofi^=N) , 

E^n^ and c^^^ m = [(:'^''^^]m where j = 1, . . . , N + I, rn. = 1, 2. .... + 1 as eigenvalues and eigen- 
vectors of H^'^^ respectively, E m cind Cm — [Cm y as eigenvalues and eigevectors of H respectively. 
In other words 

if(^)c(^)^ = Si^'c^f ) and HCrr, = EmCm. (3) 

The main aim of present work is to understand the dependence of the spectrum of H^'^^ on N . 

3 The spectrum of cut momentum and coordinate oper- 
ators 

Matrix elements of the P and Q operators in the occupation number basis read 

(n I Q I fe) = sj^5n,k-i + \J ^^^^S„,k+i , {n\P\k) = ^\J^5n,k-i + \^J'^-^Sn,k+i- (4) 

^We are considering here hamiltonians with potentials being polynomials in variable Q. 
Other types of potential functions (e.g. ^) may be analyzed as well by introducing coordinate 
representation, however numerically it is more time consuming. 
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In the Hilbert spare limited to maximum of N quanta the eigenvalues of, e.g., momentum are given 
by zeros of the determinant 



In+1 = 



-r] VT 
-Vl -V V2 
-V2 -7) 



N -ri 



, where 77 = iV2X. 



(5) 



Determinant (5) is evaluated by solving recursion relation following from the Laplace expansion. 
Making a change of variables Jn = -^In we obtain 



(N + 2)JN+2 + nJN+l- Jn =0, Jl=-V: J2 = -W +!)• 



(6) 



This recursion may be solved using the generating function method. Let us define series with 
coefficients Jn 



F{x) = J\X + J2X^ + JzX^ + JiX* + . 



It follows from (6) that F{x) satisfies 

{-x + r,){F{x) + l) = -F'{x), 
with the boundary condition F{0) = 0. The solution reads 



F{x) ^exp(^- xq ) - ^ = Yl 



'7=2) ' 



(7) 



(8) 



(9) 



where Hn{x) stands for N-th Hermite polynomial. Since F{x) is analytic at a; = 0, the expansion 



(9) is unambiguous so JN = 2-^i" Then 



In = 2- ^i^HNiX). (10) 

It is clear now that the spectrum of operator P in a cut Fock space is given exactly by zeros of 
Hermite polynomials. Therefore, denoting as the m-th eigenvalue of cut momentum P'"', we 
get 



where i/jv+i(2;^+^^) = 0, m = 1,2, . . . , N + 1. 
This result will be used several times below. 



(11) 



Calculation for coordinate operator Q is very similar. Recursion relation is slightly different but 
initial conditions also change. Those two differences cancel each other and finally we obtain the same 
result as for P. Therefore, denoting as the m-th eigenvalue of cut coordinate (3^^\ we obtain 

9m = «L^+'\ where li-jv+i(2^+'') = 0, m = 1, 2, . . . , iV + 1. (12) 

Since roots of Hermite polynomials are symmetric around 0, we consider only positive ones for which 
we introduce the following enumeration 
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0=Po <Pi <P2 - -• <Pm < -•• <Piv , N-even, (13) 

2 

Pi <P2 ■•• <P^< ■■■ < , N - odd. (14) 

2 

4 The continuum limit — scaling 

Because of the continuum limit it is particularly interesting to analyze the behavior of roots of 
Hermite polynomials when N ^ oo. It is possible to obtain the asymptotic relation ( details are in 
appendix B ) 



If one naively evaluates the limit iV — » oo for fixed m one obtains limN^ao Pm = 0. This is 

m fixed 

unacceptable, because we know that Hmjv-^oo Pm = P with p =^ and p oo. It is clear now that 

continuum 

m has to depend on N as follows 

\/2N 

m=- p + b. (17) 

TT 

A prescription, that guarantees existence of the continuum limit 

lim Pm(N) = const., (18) 

is called scaling. The dependence (17) is universal, that is for a large class of observables one obtains 
nontrivial values when iV — > oo. Substituting (17) to (16) and ordering resulting expression with 
respect to powers of N we get 

Notice that b has no influence on the result obtained in the continuum limit. Nevertheless, it is clear 
that taking 6 = | gives the best convergence. 
Moreover one can put 

m= — pViV + &+-^=, (20) 
TT ^/N 

providing another parameter that controls the convergence. Now equation (19) is modified to 
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The procedure may by continued but the optimal values of the coefficients b, c are not universal, 
i.e. if we take another observable they will be different. We can see this already in the example 
(19) where coefficient b depends on parity of N. Nevertheless the limit (18) is valid for different 
observables. 

It is interesting to deal with the problem of cardinality of the spectrum of the momentum operator. 
For all N the spectrum of cut operators consists of finite number of eigenvalues but we know that 
in the continuum limit there has to be an uncountable set of eigenvalues. How those two facts can 
be brought together? According to (15,16), for large N, there are N eigenvalues (-^ positive ones 
and -y negative ones, or ^^-^ positive ones and negative ones for N even or odd respectively) 
separated by the distance It means that the spectrum becomes denser so that it is possible 

to chose such m{N) that 

Vpefl : limN^ooPm(N) = P- (22) 

Therefore the set of all roots of all Hermite polynomials Z = U?f=i spectrum{P^-^^) is dense in R. 
However it is not equal to R because Z is countable due to the fact that there is countable amount 
of Hermite polynomials. In other words elements of Z behave similarly to rational numbers in R. It 
looks as if there was something wrong because the spectrum of operator P should be continuous. In 
order to solve this paradox we use the scaling (17). Now any real number p can be obtained in the 
continuum limit so that all elements of R are reproduced. 



5 The spectrum of the cut kinetic energy 

In order to calculate the eigenvalues of a free particle we introduce the cut parity operator 



1 
-1 







(23) 



N+i columns 

A straightforward calculation shows that [(p2)(^\ E^^'] = [(p*^')^ T,'-^^] = so that (p^)'^' as well 
as (p'^^)^, represented in an eigenbasis of E'^-', splits into two blocks. 

Let A'^ be an odd number. In this case the matrices (p^)'^' and (p^'^^)^ contain two blocks ^^^^ x -^^^ 
each. We have ^ 



a: 



(Ar)s2 



B' 



(24) 



^The dimensions of those blocks are equal in this case because for odd the rank of the 
matrix SC^) is A + 1 - even. Therefore it contains the same number of +1 and —1. 
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where 



and 



,(M) 



1 
4 

4 



4 



^(2M-3)-(2M-2) 
4 



^(2M-3)-(2M-2) 
4 

4M-3 



(25) 



,(M) 



3 

4 

4 



/2^ 
4 



^(2M-2)-(2M-l) 
4 



^(2M-2)-(2M-l) 
4 

4M-1 



(26) 



(M) 



/1-2 
4 



/1-2 
4 



^(2M-3)-(2M-2) 
4 



^(2M-3)-(2M-2) 
4 

4M-3 



(27) 



Since A''_^''' — B)^'' (the + sign corresponds to odd m) therefore 



,(M) 



,(M) 



/2-3 
4 







/2-3 

4 



y'(2M-2)-(2M-l) 



^(2M 2)-(2M-l) 
4 

4M-1 _ 2M 
4 4 



(28) 



/ 1 2\N 



^(Pm)^ where m = 1, 3, . . . , A/" for iV — odd. 



When N is even, the analogous procedure gives 



(29) 



-E^m = (^P^)m = ^(Pm)^ where m = 2, 4, . . . , iV for iV - even. 



(30) 



Now it is useful to present the eigenvalues (29,30) in the following table. For example for iV < 4, 
m < 5 
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N = N = l N = 2 N = 3 



m = 1 
m = 2 
m = 3 
m = 4 



/ = ? El = 

= ? 



El = W 



El = \{j>l? \ 



2\2 



Ei 



El = l 



El 



El = l 



(31) 



We prove that fields filled with questionmarks in (31) are equal to their neighbors on the right 



Let N be an odd number. We have already shown that matrix (^P^)^^^ contains two blocks (24). 



Now, if we increase N 



iV + 1 we obtain 



N + l 
















• 




... • 











iV + 1 










(32) 



so that the block does not change. ^ It means that in the A'' + 1 cutoff, eigenvalues from this 

block remain untouched. This block corresponds to even m, therefore we have 



Ei^'^ for m = 2,4,. 



,N + 1 N - odd. 



When A'^ is even an analogous procedure gives 

£;L^+^)=i?W for m- 
That completes the whole spectrum of {^p^)^'^^ ■ First few exemplary values are 



1,3, ...,Af N-even. 



(33) 
(34) 



0.25 0.25 


0.137 


0.137 


0.095 


0.095 


- 0.75 


0.75 


0.459 


0.459 


0.333 




1.362 


1.362 


0.892 


0.892 






2.040 


2.040 


1.400 








2.762 


2.762 










3.516 



Above numbers were obtained from a program described in next section and indeed confirm (33,34). 
According to (15,16) formulas (29,30,33,34) give 



2 2A r + 3 



N - odd m - odd 



(35) 



^The matrix A 2 becomes larger because the increment N — ► N+l produces new state 
with parity S = +1. Next increment (i.e. N + l — > N + 2) will produce new state with 
parity S = — 1 ect. 
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Si^'«Y2FT5 m-even, (36) 

E^r^^'^i^-Z^ N-even m - odd , (37) 



2 2N + 5 
~ 2 2iV + 3 



N - even m - even , (38) 



Note that (17) applied to each (35,36,37,38) separately gives the expected limit Moreover we 
see that the dependence of spectrum on TV is power like i.e. slow. 

6 Applications 

In this section we use above analytic results to verify the method introduced in [7,8]. It consists 
of numerical diagonalization of finite matrices and extrapolation of results to N ^ oo. Practically, 
when one deals with fast convergence of eigenvalues it is sufficient to stop the calculations for rela- 
tively low cutoff N ( in the case of one dimensional nonrclativistic quantum mechanics the results 
for N=50 are already very accurate ). Nevertheless a problem may occur when the convergence is 
slow ( polynomial ), or when numerical calculations are time consuming even for low N. 

One of aims of this work is better understanding the case of a free particle which has the former 
feature. The later situation occurs everytime when there are higher dimensions. Models discussed in 
[7,8] have both of those difficulties, therefore it is crucial to understand analytically the asymptotics 
of the spectrum for large N. We expect that the power-like behavior in N is characteristic not only 
for the spectrum of a free particle but also it occurs in every scattering problem because in those 
rases the asymptotics of wave functions is the same as for a free particle so that the asymptotic 
momentum may be properly defined. 

6.1 Quantum mechanics on a computer 

Let us discuss in details the implementation of the method [7,8] in the computer code. Consider 
quantum system with D degrees of freedom with D creation and annihilation operators. One can 
construct the whole orthogonal basis from the vacuum state | 0) 

I n.n.,...,n.)^ 10). (39) 

Each state in a cut Fock space, decomposed in this basis, is represented as a list in Mathematica 
program 

V 

|p,{ai,a2,...,ap},{{ni,n2,...,ni,},{ni,'M2,...,Mo},...,{<,n5,...,n^}}|. (40) 

The first element of this list specifies the number of basis vectors, that the state | V') is decomposed 
on . The second element of the list is a list of coefficients of this decomposition. Basis vectors are 
represented in the third element of this list. For example 

a I 0, 1) + fej 1, 0) + c I 1, 1) ^ {3, {a, 6, c}, {{0, 1}, {1, 0}, {1, 1}} }. 
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The creation and annihilation operators 

ak- ttk \ ni,n2,. . . ,nk,- ■ ■ ,nD) = Vnk \ ni,n2,. . . ,nk - 1,. ■ ■ ,nD), (41) 

Sfe ■ all ni,n2, ■ ■ ■ ,nk, ■ ■ ■ ,nD) = Vnjc + 1 | ni, n2, . . . , + 1, . . . ^no)- (42) 
have the following action in the list representation 

flfc I V') — > [p, {\l nlai, \lr^),an}, {{nj, . . . , - 1, . . . , rir,}, {n?, . . . , - 1, . . . , n|,}, . . . , 

...,{<,...,<-!,.. .,n^}}}, (43) 

and 

I V*} — * [p, {y + lai, ■ ■ ■ , + lap}, {{n\,..., ni+l, . . . , nL}, {n?, . . . , n^+l, . . . , nl,}, . . . , 

...{<,...,< + !,. ..,n^}}}. (44) 

In order to evaluate the matrix representation of any observble we define procedures which add and 
multiply on arbitrary state by a complex number as well as scalar multiply states. For example 

IV'>^{2,{1,2},{{0,0},{0, 1}}}, 
10) ^{2, {1,1}, {{0,2}, {0,1}}}, 

then 

I V'>+ 1 4>) {3, {1, 3, 1}, {{0, 0}, {0, 1}, {0, 2}}}, 
2 10) ^{2, {2, 2}, {{0,2}, {0,1}}}, 

and 

(V I 0) ^ 2. 

Adding lists is simply adding those coefficients of the decomposition (40), that have the same basis 
vectors. If decompositions of | V) and | ip) have different basis vectors then the sublist consisting of 
basis vectors has to be extended accordingly. 

The procedure of multiplying the state by a number reduces to multiplying the list of coefficients by 
this number. 

Scalar multiplication {ip \ (j>) reduces to a search for common basis vectors occurring in decomposition 
of I tp) and I 0) . Afterwards proper coefficients and their complex conjugations have to be multiplied. 

These rules allow to automatically represent any operator in a cut basis (39). 
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6.2 Numerical diagonalization 



Here we compare numerical data and analytic results of section 3for 

al) eigenvalues of p''^' evaluated by the program described in 6.1. ( according to section 3 they are 
exactly the roots of Hcrmite polynomials ), 
a2) the asymptotic form (16), 



0.1 * 
. 6 
0.5 * 



♦ al 

★ a2 



0.4 

I 2 I V%*hO ^ 



P m=2 



♦ al 

★ a2 



0.8f 



P\^3 

3.5 y 

3 

2.5 a* 



1.5 



♦ al 

★ a2 



-171$|4T.0 



Figure 1: Asymptotic, a2), and exact, al), behavior oipm^ for m=l,2,3. 



Figure 1 presents the comparison of cases al and a2 for m=l,2,3. The approximate value is obtained 
from (16) by taking only the leading term 



We see that there is a good agreement between exact and approximate values even for low A'^, and 
it gets worse for higher m where next terms of the expansion of (16) are important. 



6.3 Continuum limit on a computer 

Here we want to obtain dispersion relation that is the dependence of the energy on momentum E{p). 

2 

Obviously we know that E{p) — ^ but it is only because we are able to solve Schrodinger equation 
for a free particle. However one has to put himself in a situation where there is a certain set of 
eigenvalues and no information about the dispersion relation is available. In other words the 
question is how to obtain unknown a priori function E{p) by means of eigenvalues E^? In order to 
do this one has to make m dependent on N : m = m{N,p) such that the limit 

lim =£;(p), (46) 

AT— »oo 

is not trivial that is E{p) < oo and E{p) / 0. Note that (46) automatically requires the set 
{E^ : m, e N} to be dense in E{R) . In case of a free particle (i5(R) = [0, oo)) we can even 
construct this set (squares of roots of Hermite polynomials) however it is a general property of any 
operator with continuous spectrum. This is exactly the reason why E^ depend on N as a power 
rather then exponentially. 

Let us emphasize that we don't have to know the dependence E{p) to evaluate m{N). This 
is because the relation m{N) was established on grounds of the condition that there has to exist 
the continuum limit for the momentum, so that any other operator commuting with P will have 
the same scaling. We will analyze in details the case of a free particle in nonrelativistic quantum 
mechanics but another example may be Dirac equation where we expect that the scaling law (18) 
will give E{p) = y'M^ + pP. Therefore the scaling in (46) has to be the same as for momentum 
operator, that is 

, VM 1 

m{N,p) = P+-. (47) 

TT Z 
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However in formula (47) we have to introduce a certain change 



m{N,p) — ^ 2m{N,p) ^2- p+1, (48) 

TV 

because the scaling (48) is meant for positive eigenvalues of operator P^^^ only. Let us consider an 
example of = 7. The spectrum of operator P'*' consists of roots of Hs{x), so that we have 8 roots 
where 4 of them are positive and 4 are negative. 



* * * * 9 • • • • — »- 

Now, if we square them the spectrum becomes positive and the numeration of eigenvalues changes 
as follows. ^ 



9 * • * • * • * 

For example, the eigenvalue that we used to number as the first one will now have the index m = 2, 
the eigenvalue that we used to number as the second one will now have the index m = 4 ect. There- 
fore the formula (47) has to be rescaled as in (48). 

According to (48), eigenvalues are analyzed by fixing any momentum value p and writing 

down the value ^-miTv^iai p) '^'^'I'S'"^ Nmax is the highest N in computer calculations. Then we change 
the momentum value and repeat the procedure. In this way one obtains an approximate ( because 

2 

of limited value oi N ) dependence E{p), which should reproduce ^ for a free particle. However the 
problem concerning the formula m{N,p) occurs because m is not a natural number. We circumvent 
this by taking an integer part ( INT ) of Eq.(48), so that the matrix index is INT{m{Nmax,p)), 
where is an even number. The convergence of those elements was checked in Mathematica for 
p= 1,2, 3,..., 37 (e.g. Figure 2) 

This behavior can be understood as follows. If one plots the dependence of on N ( p is fixed 
), one obtains (35-38) a hyperbola. The lower index m = 1 specifies the first eigenvalue. The upper 
index enumerates the cutofi^. If we plot the dependence of E2 on N, we get another hyperbola ect. 
Finally the plot of E^ is a set of hyperbolas on a plane ( see Figure 3). The scaling that we have used 
previously means that from each hyperbola we are taking only one point in such way that in the limit 
of large N a constant value is reproduced. Why on those figures we see cut hyperbolas instead of 
points? This is because we had to introduce the INT procedure which is equivocal. In a consequence 
it is possible that for different cutoffs ( say N and N' ) there is INT{m{N,p)) = INT{m{N' ,p)). It 
means that points {N,INT{m{N,p))) and {N' , INT(rn{N' ,p))) are on the same hyperbola. Even- 
tually N will be large enough so that the INT operation notices the difference and the point "jumps" 
to next hyperbola. Let us also note that the scaling (18) is an asymptotic law hance for low N the 
behavior of E^^j^ ^-^ may vary for different values of p. This effect accounts for the different behavior 

■^Since roots of Hn(x) are symmetric around the origin, their squares will give double 
degeneracy. Hance for the free particle dots and stars should be on the same point however 
in general it is not the case. 
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in Figure 2a and 2b. 

The dispersion relation extracted in this way is presented in Figure 4. This result has no error 
because all eigenvalues are precisely evaluated hance any statistical interpretation is meaningless. 
The tangent coefficient 0.5365 is different from 0.5 but we did expect that because it is a numerical 
result obtained on grounds of limited cutoff. Moreover we had to introduce the INT operation. In 
a consequence we had to choose only one point from cut hyperbolas. It is a source of a new error 
which gets smaller while the cutoff increases. 

Therefore Figure 4 confirms that we can obtain the dispersion relation from the knowledge of the 
spectrum of a cut hamiltonian. 
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4 6 8 10 

Figure 4: Reproduced dispersion relation. 



7 Bound states versus scattering states 

In this section we stress the difference between localized and nonlocalized states. It follows from 
simple algebra ( see Appendix A) that ^ 

which means that the spectrum of cut operators converges towards the spectrum of operators in 
infinite Hilbert space. Moreover one can tell how fast is the convergence because from (49) it is 
clear that the convergence > Em is governed by the behavior of the at large j. Note 

N »oo 

that in (49) arc the exact components of eigenvectors of H. This is exactly the result we were 
anticipating because the difference between localized and non- localized states lies in components c^. 
Therefore one can numerically judge weather the state is bound or not on grounds of the behavior 
of the eigenvalues of cut operators only. 
For the case of a free particle one can obtain exactly 

cj = (n|fe)= / dx{n\x){x\k) = [ dxij" ° {x)e''"' , (50) 
Jr J r 

where H.O. stands for harmonic oscillator 

°(x) = -=i=/?„(x)e-=^'/^ (51) 
V2"n!7r 

Integral (50) is evaluated with the aid of some analytic properties of Hermite polynomials, what is 
^The notation is explained in section 2 
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presented in Appendix C. The result is 



Figure 5 is an example of (52) ior E = = 1000. 



(52) 



2-10 



1-10 



242 



-1-10 



-2-10 



242 





fl 






^0 o 











00^ 



Figure 5: Components of the eigenvector (E=1000) for free particle. 



Asymptotic behavior of the envelope is ( see Appendix C ) | |« y ^ which is indeed power like . 
Similar calculations for discrete spectrum are not known, so one is left with numerical data instead. 
Figure 6 presents components of eigenvector corresponding to the first (the lowest) eigenvalue of 
anharmonic oscillator, as well as the convergence of the first eigenvalue. 

In this case the behavior of is completely different from one shown in Figure 5. One sees that 
varies in the same ( exponential ) way as c^. In other words, the behavior of eigenvalues E^ 
with the cutoff N distinguishes wether the state is bound or not. 



8 Conclusions 

The main purpose of this paper was to prove that the method proposed in [7,8] enables to distinguish 
numerically weather the state is localized or not. This aim and related problems have already been 
investigated [9-13]. This distinction is an important issue while studying supersymmetric models 
(D=10 SYMQM) where bound states exist among dense number of scattering ones [1]. Therefore 
one has to reanalyze quantum systems from the very beginning in a new manner. Starting from the 

•^Eq. (52) can be obtained independently in a shorter way. Notice that is a Fourier 
transform of ip^'^'{x) which is the solution for hamiltonian H = Ip^ + . The Fourier 
transformation switches x with p but H is symmetric in those variables so the Shrodinger 
equation in momentum representation is the same as in coordinate representation. Therefore 
the solution for harmonic oscillator in momentum representation is of the same form (up to 
a multiplication factor) as the solution for harmonic oscillator in coordinate representation. 
The connection between those two solutions is given by Fourier transform hence coefficients 

are of a form (52) . 
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Figure 6: Components of the eigenvector (c^^i) and the the convergence of 
eigenvalues (E^^i) for anharmonic oscillator. 



calculation of spectrum of cut operators Q ' ',P^ 'one realizes that eigenvalues of those operators 
are exactly equal to the roots of Hermite polynomials. Next, we conclude that in order to recover the 
continuum limit one has to introduce scaling m{N). The validity of the scaling law in the hamiltonian 
of a free particle was rigourously proven in section 5 and numerically tested in section 6. As a result 
one reproduces the dispersion relation from an information about a spectrum of a cut hamiltonian. 
It is expected that the same scaling may be applied for a set of hamiltonians commuting with P or 
under weaker assumptions, namely those for which P can be defined asymptotically. The scaling 
in higher dimensions is important because of the occurrence of scattering states (e.g SYMQM D=2 
systems). The formula (18) is expected to be valid in those cases because they are described by 
quantum mechanics of a free particle in color dimensions. In this case the coefficient in (17) may be 
different, however (18) is claimed to be applicable all the time. In particular D=2, SU(2) SYMQM 
[10] is free and it has been found [14] that the system requires (17) exactly to recover the continuum 
limit. Recently a new possibility to speed up the numerical approach in D=4 has occurred Jll]. The 
naive diagonalization of the Hamiltonian in the whole cut Hilbcrt space was abandoned and replaced 
by the language of rotational invariance. The new approach can be extended to higher dimensions 
as well. 
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10 Appendix A 



Here we derive the formula (49). Let us start with eigen equation Hcm = EmCm. where H is an 
operator and Cm its eigenvector. Writing it in the matrix form 



(53) 





hi JV+1 




r 1 

*--m 




r 1 




flN JV+1 






~ Em 




flN+l 1 • • • tlN+l N 


flN+l N+1 











and rewriting for first N components only one obtains 

,^1 n\ N+iCm 



H 



(N) 



r 1 






+ 

















r 1 


— Em 













Now complex conjugate (54) and multiply it by c^^'n from the right side 



(54) 



[ ^ m 



I.JV -1 jj(N) 



N+i 

m ? • • • : 



j=l ft N,N+iC m J 



so that 



p r *i *Af 1 



JV JV oo JV 

piN)Sr^ *j (N)i Sr V h*-r*^-F Vr*^' 
i=l i=l 3=iV+l <=1 



1^ 



(55) 



(56) 



JV JV oo 



[Em — En ') C '„ — /l ijC n, 

i=l 1=1 j=JV+l 



(57) 



it is non trivial to realize that above equation means that YliLi ^*rnC^'^^n — * ^mn thus one can omit 
the n index and write 



EN oo 7 



c*W' ^^^^ 

Of course this derivation is for the case with discrete spectrum (discrete index m ) nevertheless for 
continuous spectrum the same calculations give 



<=1 i=l j=N+l 



J 



(59) 
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where 

^(iv)c(iv)^ = ^W^W HCE = ECE EGR. (60) 

This case is discussed in details in sections 6 and 7. 

11 Appendix B 

In this appendix we derive the asymptotic form of the zeros gJJ:, of the Hermite polynomial Hn{z). 
When n is an even number they may be obtained using the following relation [14] 

H„iz) = {-lf2"^n\L-^^z^), (61) 

where n is an even number, L"(«^) are the generalized Laguerre polynomials with parameter a (in 

out case a = — ^). Let ,tm,a and denote the m-th positive root oi Hn{z), Ln°'{z) and Ja(«). 
One has [14] 

a,„ = ^ (l + ^^^^1^2^) + 0{n-% (62) 
where ^ 

fcf ,a = 2 2~ iAi,af = tm,c,- (63) 

For a = — I we obtain J_i{z) — ■^J^cos{z), therefore ja,m = n{m — i) where m = 1, 2 ^ 
and kn^a = ^ + j so 



/ n \2 



2n + 1 \ 3(2n + 1) 

Let us define ( m is fixed ) 



^("^•= 2n+l ^+ 3(2n+l)2 = n + + + ' ' ' ' ^^^^ 



we have 



V/(n) + 0(n-5) = v7W^l + ^|^ « ^/(n)(l + io(n-*)) = ^/(n)(l + 0(n-*)), 

(66) 



C = VTR + y7RO(n-^) = VTR + 0{n-'-'), (67) 

finally 



Tim V + 3(2n+i)^ +^(" ^^^^ 
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When n is an odd number there are [10] analogous relations 



and 



,2 jm,a:'' I ^ 2(q^ 1) + jm,a:^ \ Q/ -5-. 



4k 
where 

fe 't-l ^ = ^ I ^ and (^m.ct) ~ im^a • 

In this case Ja{z) = Ji (z) = •\/^sin(«) so ja,m. = wm where m = 1,2...,..., and 

n I 1 
2 ' 4- 



§ + 1. Analogous calculations give 



^"--^2^^/ 3(2n + l)2+^^" ^• 

12 Appendix C 

Here we evaluate the integral 

J R 

It follows from three properties of Hermite polynomials [10] that 

Hn{x + y) = ^Y1 l2,]Hm{V2x)Hn-m{V2y), 

m=0 V / 

J7n(a;) = ^ / dt{x + itTe-* , 
/ da;i?n(a;)-ffm(a;)e""'^ = 2'^n\y/T^5nm- 

J R 

After substituting (75) to (73) and changing the variables x — > a; + ifc we get 

/4fe) = -^2^e-'='/V| dte-'" j dxe-^" {i±^+ix 
Using (75) once again we obtain 

74fc) = e-'='/^2^i" j dte-*"H^ 
Finally substituting (74) to (78) and using (76) we get 



_k 



I^{k) = e-'='/'V2r J2(j H^-m{k)2"'m\V^5mO = e-'''/^V2^i"H„{k), 
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therefore 



It is straightforward now to estimate components 




(80) 



1 



In{k) \< 



1 



/n(0) |< 



1 



V2^ 1 H„iO) I . 



(81) 



Ce 1 = 



Since -ff2n+i(0) = hance j 4"+^ |= 0. On the other hand ii'2n(0) = (-l)"-^^^ therefore j c|" ]< 
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. Finaly according to Stirling formula one obtains 




(82) 
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